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We present the results of molecular dynamics computer simulations 
of a simple glass former close to an interface between the liquid and 
the frozen amorphous phase of the same material. By investigat- 
ing Fs{ci,z,t), the incoherent intermediate scattering function for 
particles that have a distance z from the wall, we show that the re- 
laxation dynamics of the particles close to the wall is much slower 
than the one for particles far away from the wall. For small z the 
typical relaxation time for Fs(q, z, t) increases like exp(A/ {z — Zp)), 
where A and constants. We use the location of the crossover 

from this law to the bulk behavior to define a first length scale 
z. A different length scale is defined by considering the Ansatz 
F,(q,z,t) = F^''^(q,t) + a(t) exp [-(^/^(t))'^^] , where a(t), ^W, 
and (3{t) are fit parameters. We show that this Ansatz gives a very 
good description of the data for all times and all values of z. The 
length ^(t) increases for short and intermediate times and decreases 
again on the time scale of the a-relaxation of the system. The max- 
imum value of ^{t) can thus be defined as a new length scale ^max- 
We find that z as well as ^max increase with decreasing temperature. 
The temperature dependence of this increase is compatible with a 
divergence of the length scale at the Kauzmann temperature of the 
bulk system. 



1 Introduction 

The idea that the slow dynamics of supercooled liquids is related in some way to the 
existence of domains in which the djTiamics of the particles is cooperative is an old one 
(Adam and Gibbs 1958). The size of such domains is supposed to grow with decreasing 
temperature and can in turn be used to explain the slowing down of the dynamics. For 
a long time it was not clear at all whether or not the idea of such domains has any 
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counterpart in real glass formers or whether it is just an appealing theoretical concept. 
Only in recent years a variety of novel experimental techniques yielded results that gave 
evidence that such domains do indeed exist in the supercooled state (e.g. Cicerone et 
al. 1995, Ediger 2000, Hempel et al. 2000, Leheny et al. 1996, Richert 1996, SiUescu 
1999, Tracht et al. 1998, Yamamuro et al. 1998). Also computer simulations of two 
and three dimensional model systems give evidence that in the supercooled state the 
relaxation dynamics of the particles is very cooperative (Dohwa and Heuer 2000, Donati 
et al. 1999a/b, Perera and Harrowell 1999, Yamamoto and Onuki 1998). Progress has also 
been made on the theoretical level in that it was shown that such domains are related to 
the nonlinear susceptibility of a certain four point correlation function (Franz and Parisi 
1998, Franz et al. 1999, Franz and Parisi 2000). 

In the present paper we use the fact that the relaxation dynamics of a particle in a 
confined system depends strongly on its distance z from the wall (Scheidler et al. 2000, 
Scheidler 2001). By investigating how this infiuence depends on z and time t we thus can 
determine a dynamical length scale and study its dependence on temperature. 

2 Model and Details of the Simulations 

The system we study is a binary (80:20) mixture of Lennard- Jones particles, i.e. with an 
interaction potential Vaisir) — ^^aisli'^ap/fY'^ ~ ('^a/?/'")^]- This potential was truncated 
and shifted at 2.5(Tq;/j. The interaction parameters are uaa = 1-0, eAA = 1-0, ctab = 0.8, 
eAB ~ 1-5, cbb = 0.88 and eBB = 0.5, where A and B label the type of particle. In the 
following we will discuss the results in terms of reduced units, using ctaa and caa as the 
units of length and energy (setting the Boltzmann constant ks = 1), and {ma\j^/ AScaaY^^ 
as the unit of time. (Here m is the mass of the particles.) 

In the past the static and dynamical properties of this system in the bulk have been 
analyzed carefully (Kob and Andersen 1995a/b, Gleim et al. 1998) and it was found 
that its dynamics is described very well by means of the so-called mode-coupling theory 
(Gotze 1999) with a critical temperature around Tc — 0.435. Thus this is the temperature 
at which the dynamics of the bulk system slows down significantly. A further relevant 
temperature for this system is the Kauzmann temperature Tk which has been estimated 
to be around 0.3 (Coluzzi et al. 2000a/b, Sciortino et al. 1999). 

The confined liquid studied in the present work was set up in the following way: We 
first equilibrated a system of 3000 particles in a box of dimensions = Ly = 12.8 
and Lz — 15 by applying periodic boundary conditions in all three directions. After 
equilibration we joined a copy of the configuration in the positive z-direction and one in the 
negative ^-direction. The particles of these two copies that were in the layer (perpendicular 
to the 2;-axis) of thickness 2.5 closest to the original system were subsequently frozen and 
constituted the wall for the confined system (2.5 is the range of the interaction). In 
this way we thus generated a sandwich geometry in which a fiuid system of extension 
Lz was confined by two walls of dimension x Ly x 2.5. (In the direction of 
the X and y axis we kept the periodic boundary conditions.) The width of the film 



2 



was sufficient to make sure that the dynamical bulk properties are realized in the center. 
Apart from the two mentioned walls we also added a hard-core potential at 2; = and 
z = Lz to avoid a penetration of particles into the wall. It is now straightforward to show 
that with these types of boundary conditions the equilibrium structure of the confined 
liquid is identical to the one of the bulk liquid. Hence it is not necessary to equilibrate the 
confined system after the introduction of the wall. (Note that this procedure was done 
for all temperatures investigated and hence the structure of the wall depends slightly 
on temperature.) Moreover the energy of a particle near to the wall is statistically the 
same to the one of a particle far from the wall: in general the presence of such a wall 
cannot be discovered by looking to quantity measured at one time, but influences only 
the correlation of quantities measured at two different times. 

The lowest temperature investigated was T = 0.5 for which production runs of 16 
million time steps were used. Note that the length of these production runs is more than 
ten times longer than the one needed in a simulation for the bulk (Kob and Andersen 
1995a/b), since the relaxation of the system close to the wall is so much slower. Finally 
we mention that in order to improve the statistics of the results we averaged them over 
8-16 independent runs. 



3 Results 

The dynamical quantity that we investigate is the incoherent intermediate scattering func- 
tion Fs(q, t), i.e. a space-time correlation function that is important from the theoretical 
as well as experimental point of view (Hansen and McDonald 1986). Due to the presence 
of the walls the system is no longer homogeneous or isotropic and hence we have to con- 
sider the 2;-dependence of Fs(q, t). For this we investigate the following generalization of 
the incoherent intermediate scattering function: 

Fs{ci,z,t) = — (^X:exp(zq- (r,(t) - r,(0)))5(z,(0) - z)j (1) 

where the 5-function on the right hand side makes sure that one considers only particles 
that at time t = were a distance z away from the wall. The wave-vectors q we consider 
in the following are parallel to the walls and their x and y components are chosen such 
that they are compatible with the periodic boundary conditions. In Fig. |T] we show the 
time dependence of Fs{c[,z,t) for various values of z. From this figure we see that the 
relaxation dynamics of the particles close to the walls (curves for ^ 0) is much slower 
than that for particles in the middle of the film (curves for z ~ 7.5). In an earlier paper 
we have shown that the typical relaxation time for particles that have a distance z from 
the wall, r{z,T), increases like 

t{z,T) ocexp [A/ {z - Zp)], (2) 

where A 7.5 and Zp ^ —0.5 are constants and the proportionality factor depends on T 
(Scheidler et al. 2000, Scheidler 2001). This law is valid for small values of z, i.e. close to 
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the wall. For large z we recover instead the bulk relaxation time T\^^x\^J(^\ The distance 
from the wall at which one finds the crossover from the law given by equation (0) to the 
constant bulk value rbuik(7') can be used to define a length scale z (Scheidler 2001). Thus 
z is the length scale over which the wall infiuences the dynamics of the particles and below 
we will discuss how this scale depends on temperature. 

A second approach to investigate the spatial dependence of the structural relaxation is 
to consider the whole t- and z-dependence of the intermediate scattering function within 
one Ansatz. From figure H we see that the z-dependence of F<j(q, is very smooth. 
In particular it is clear that for z ^ \ this function has to converge to the intermediate 
scattering function for the bulk, F]'"^'^(q, t\ Therefore it is reasonable to make an Ansatz 
of the form 

F,(q,^,t) = F^i'^(q,t) +a(t)exp [-(^/^(t))^^] , (3) 

i.e. that the whole z-dependence of Fs(q, is given by a stretched exponential. Here 
a(t), ^(t), and /3(t) are unknown functions of time that can be found from fitting at a 
given time the z-dependence of -F<,(q, z, t) with the functional form given by equation (^. 
Note that F]'"^'^(q, t) is known from the simulations of the bulk and is hence not a fit 
parameter. 

That the Ansatz (|^) is indeed able to describe the z-dependence of Fs(q, very 
well is demonstrated in figure ^ where we show F5(q, 2;, t) — F]^""^(q, t) as a function of 
z for times that span the range from microscopic times to the time of the a-relaxation, 
i.e. of the final decay of the correlation functions. We also mention that other simple 
functional forms, e.g. if the stretching parameter /5 is set to 1.0, do not give satisfactory 
fits (Scheidler 2001). Nevertheless such a purely exponential Ansatz leads to a very similar 
dynamical length scale. 

In figure |^ we show the time dependence of the stretching parameter /9(t) for all 
temperatures investigated. From this figure we recognize that at low temperatures /?(t) 
shows at short times (~ 3 time units) a rapid decrease from values around 1.5 to a value 
around 1.2. Subsequently it remains at this latter value for a time span that corresponds 
roughly to the /3-relaxation, i.e. the time window during which the correlation functions 
for small z remain close to the plateau and which extends at the lowest temperature over 
roughly three decades in time (see figure |I]). For longer times /3(t) increases again to values 
above 1.4. For intermediate and high temperatures the time dependence of /3(t) remains 
qualitatively the same as the one we just described. The main difference is that the 
length of the plateau at intermediate times is shortened significantly, in agreement with 
the fact that also the correlation functions show only a short plateau at high temperatures 
(Scheidler 2001). Note that the fact that at intermediate times /?(t) is larger than 1.0 
shows that the ^-dependence of the Fs(q, 2;,t) is rather strong, i.e. it is not just an 
exponential dependence. Additionally we observed in thin films that the infiuence of two 
walls on the dynamics of the particles between them is not just the superposition of two 
laws of the form given by equation (^. This is thus evidence that the slowing down of 
the dynamics is related to the presence of a strongly non-linear process such as, e.g., the 
non-linear feedback process of mode-coupling theory. 
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In figure ^ we show the time dependence of the parameter ^ [t) for all temperatures 
investigated. We see that for short and intermediate times this length scale increases, 
attains a maximum around a time that corresponds to the a-relaxation time of the bulk 
system, and then starts to decrease again. The value of ^ at its maximum can thus be 
used to define a new length scale C,max{T). Note that the time at which this maximum 
is attained is roughly in the time window in which (3{t) has a minimum, i.e. is relatively 
small. This implies that at this time the infiuence of the wall on the dynamics of the 
particles in the fiuid extends over the largest range. 

Finally we discuss the temperature dependence of the length scales z and ^max- These 
quantities are shown in figure |^ as a function of temperature. We see that both length 
scales increase with decreasing temperature but that this increase is relatively modest. In 
the temperature interval 0.5 < T < 1.0, where the relaxation times of the system in the 
bulk increases by about a factor of 500, this increase is only around 2.5 for the case of ^max 
and a factor of 4 for z. Due to the slow relaxation of the particles close to the wall it is 
presently not possible to study the relaxation dynamics of these particles at significantly 
lower temperatures and hence it is also not possible to determine neither z nor .^max at 
low temperatures. Due to this relatively weak increase of the length scales it is difficult 
to make precise statements on whether or not these length scales diverge or not and if 
yes, at which temperature this would be expected to happen. Furthermore it is presently 
not even clear whether such a divergence occurs at the same temperature, since the ratio 
5/^max depends on temperature (see inset of figure |^). Since previous simulations have 
identified two relevant temperatures, the critical temperature Tc of mode-coupling theory 
at Tc = 0.435 and the Kauzmann temperature Tk at around Tk = 0.29 (Coluzzi et al. 
2000a/b, Sciortino et al. 1999) we have tried to see whether the increasing length scale 
is compatible with a power-law divergence at one of these two temperatures. Whereas it 
seems that no such divergence occurs at Tc, the data is indeed compatible with a critical 
behavior at Tk with an exponent of the power-law around 1, both for z and ^max- This 
result is in agreement with a theoretical approach to describe the temperature dependence 
of the size of cooperatively rearranging region in the concept of configurational entropy 
(Huth et al. 2000) . However, it is also possible that the divergence occurs only at T = 
since also this temperature is compatible with our data. Hence we conclude that although 
we are able to identify a growing length scale it is presently not yet possible to give a 
definite answer regarding the precise temperature dependence of the scale. Therefore 
simulations at lower temperatures as well as some theoretical guidance on this matter 
would be highly desirable. 
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Figure 1: Time dependence of the generalized incoherent intermediate scattering function 
from equation (1) for the A particles for T — 0.55. The curves are for q — 7.18, the 
location of the main peak in the static structure factor. The different curves correspond 
to different distances z from the wall (see label of curves) . 
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Figure 2: The time dependence of Fs(q, z, t) — F^"^^{(\, t) is shown to check the vahdity of 
the Ansatz given by equation (3) (the value of q is the same as in figure 1). The symbols 
are the data from the simulations and the smooth curves are the fits with the functional 
form given by equation (3). 
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Figure 3: Time dependence of the stretching parameter (3{t) from equation (3) for all 
temperatures investigated. 




Figure 4: Time dependence of the scale ^{t) from equation (3) for all temperatures in- 
vestigated. The value of the maximum of the curves is used to define the length scale 

max- 
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